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ABSTRACT 

We study the evolution of QCD phase transition-generated magnetic fields in freely decaying MHD 
turbulence of the expanding Universe. We consider a magnetic field generation model that starts 
from basic non-perturbative QCD theory and predicts stochastic magnetic fields with an amplitude 
of the order of 0.02 fiG and small magnetic helicity. We employ direct numerical simulations to model 
the MHD turbulence decay and identify two different regimes: "weakly helical" turbulence regime, 
when magnetic helicity increases during decay, and "fully helical" turbulence, when maximal magnetic 
helicity is reached and an inverse cascade develops. The results of our analysis show that in the most 
optimistic scenario the magnetic correlation length in the comoving frame can reach lOkpc with the 
amplitude of the effective magnetic field being 0.007 nC. We demonstrate that the considered model 
of magneto-genesis can provide the seed magnetic field for galaxies and clusters. 
Subject headings: primordial magnetic fields; early Universe 



1. INTRODUCTION 



The origin of the observed 
m galaxies and clusters of 



magnetic fields (MFs) 
10- ^ - 10"^ Gauss 
(G) remains a matter of debate (|Beck et al.l 119961 : 

Widr ow 2002; Vallee R ecently several differ- 

ent groups (Neron ov fc Vo vk 2010':'Tav ecchio et al.ll201C ; 



[Polag et all .201 It iTavlor et al 2011; .Huan et al.ll2011 ) 
reported the detection of a lower bound on a large- 
scale correlated MF amplitude of the order of 10 — 
10~^^ G, or possibly two orders of magnitude smaller 
(jDermer et al.ll2011l :lTakahash i et al] [2Q12) at Mpc scales 
through blazar observations. One of the possible ex- 
planations of the large-scale correlated MF assumes the 
presence of a seed primordial magnetic field (PMF) which 
was generated during or prior to the radiation dom- 
inated epoch. This MF should satisfy several condi- 
tions: (i) The PMF should preserve approximate spa- 
tial isotropy, it has to be weak enough when its energy 
density can be treated as a first order of perturbation; 
(ii) The PMF should be smaller than the MF in galax- 
ies by a few orders of magnitude at least, since dur- 
ing structure formation PMFs get amplified; (iii) Since 
the PMF energy density pB contributes to the radia- 
tion field, the big bang nucleosynth esis (BBN) bound im- 
plies Q^Bhl = pb/Pct < 2.4 x 10"^ rasso fc RubinsteinI 
I2OOII ). where pcr is the critical density,-^ and ho is the 
Hubble constant in units of 100 km s~^ Mpc~^. 

The possible origin of the PMF from the two major 
cosmological phase transitions, the electroweak phase 
transition (EA VPT) and the QCD phase transition 
(QCDPT) fsee lGrasso fc R.ubinsteinl[200TI : IWidrQwll2002l : 



^ The ratio of ps to the energy density of the radiation prad is 
constant during cosmological evolution if the PMF is not damped 
by a MHD (or other) process and therefore stays frozen into the 
plasma. 



I Widrow et al.l [20121 for reviews) is of particular impor- 
tance for cosmology. Because of the larger scale of the 
resulting seed magnetic field and the nature of the QCD 
bubble walls during a first order QCDPT, it is more likely 
that the QCDPT rather than the EWPT produces a 
PMF that accounts for the observed galactic and clus- 
ter MFs. 

In this paper we consider one of several possible mecha- 
nisms of PMF gener ation. In p articular we re-address the 
model proposed bv iiKisslingerl (j2003), in which the PMF 
is generated via QCD bubble collisions. We consider the 
coupling of this initial PMF with the QCD plasma, and 
study the dynamics during the expansion of the Universe. 
The main parameters of the described model are given by 
the QCDPT temperature = 0.15 GeV and the number 
of relativistic degrees of freedom = 15. The interac- 
tions between the PMF and the QCD plasma is stud- 
ied through numerical MHD simulations using the Pen- 
cil Code (see http:/ /pencil-code.googlecode.com/). We 
discuss observational signatures of such a QCDPT PMF, 
including observed MFs in galaxies and clusters. We em- 
ploy natural units with h = 1 = c and gaussian units for 
the MHD formulation. 

The outline of the paper is as follows: In Section 2 we 
describe the PMF generation model. In Section 3 we de- 
termine the spatial and temporal characteristics of the 
generated PMF. The results of our analysis, including 
the dynamics of the PMF, are presented in Section 4, 
where we discuss the resulting MF in galaxies and clus- 
ters. Conclusions are presented in Section 5. 

2. MAGNETIC FIELD GENERATION MODEL 

In contrast to the EWPT, the QCDPT involves the 
treatment of QCD, which, unlike the electroweak theory, 
is non-perturbative. Therefore a valid theory starting 
from basic QCD theory, rather than a model, must be 
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able to treat non-perturbative QCD. In lKisslingerl (|2QQ3f ), 
instantons form gluonic bubble walls and it is the interior 
gluonic wall that leads to the magnetic seed described 
below satisfy that criterion. In this early work the main 
interest was the prediction of polarization correlations 
in cosmic microwave background radiation (CMBR). As 
we can see below the magnitude of the resulting MF is 
too small for current CMBR observations, but it might 
be measured in the future. Because this scenario starts 
from basic non-perturbative QCD theory and success- 
fully predicts a primary magnetic field which has the 
overall properties that are promising for the PMF, we 
use it in our present work. 

In this section we briefiy describe the PMF sce- 
nario proposed by 'Kissl ingerl ()2QQ3f ). In Section 12.11 
we review the magnetic field and helicity density 
created during the QCDPT, and in Section 12.21 we 
give values of these quantities at the present time. 
Recent lattice QCD studies have shown that the 
QCDPT i s first order , so b ubbles form and col- 
hde fsee [Pe Risi et al.l [2QQ8I: iFraga fc Misherl [2QQ1 



Bessa et a l.'"2QQ9: "Schwarz & Stuke 2009; Schettler et alj 
20TTI: [Bpeckel & Schaffner-Bielich 2011, and refer- 



ences the rein) ; additional references are given by 
iKahniashvili et all (|2010b[ ). The first order QCDPT can 
result in the generation of a MF through two (or more) 
bubble collisions. 

2.1. Magnetic Field and Helicity During the QCDPT 

The QCD phase transition critical temperature is de- 
fined as 0.15 GeV. A gluonic wall is created as two 
bubbles collide, and a magnetic wall is formed by the 
interaction of the nucleons with the gluonic wall. The 
electromagnetic interaction Lagrangian is 



(1) 



where ^ is the nucleon field operator, A^^ is the elec- 
tromagnetic 4-potential, and 7^ are the Dirac matrices. 
In [Kisslingerl (2003) it was shown that the interaction of 
the quarks in the nucleons with the gluonic wall align 
the nucleon magnetic dipole moments, producing a B- 
field orthogonal to the gluonic wall. 

Using an instanton model for the gluonic wall oriented 
in the x-y direction (say), one obtains for Bz{x) at the 
time of the QCDPT, with T = T^, 

5,(^)=M^^^)e-^^(-^+^^)e-^^\ (2) 

where = dn — a few km =horizon size at the end of 
the QCDPT {t and M'^ = 0.2 fm. 5?^^, the 

magnitud e of the M F within the wall of thickness is 
see Kissl ingeiil200l 



Si 



(QCD) 



CAqcd 2M„ 



X {*(J2l75*) , 



(3) 



where Aqcd — 0.15 GeV is the QCD momentum scale, 
75 = Z7^7"'"7^7^, and (721 = ^7271 = ^7^7"*". A simi- 
lar form h ad been derived earli e r usin g; the domain wall 
model of iForbes fc Zhitnitskvl (|2000[ ). The value for 

^(QCD) fo^^^ 



^(QCD) 



0.39-A^CD - 1-5 X 10-^ GeV^ 

TT 

2.2 X 10^^ G . 



(4) 



The asterisk indicates that we refer to the initial value 
of the MF at the time of the QCDPT. 

We now discuss the magnetic helicity crea ted during 
the QCDPT using the scenario proposed by iKisslingerl 
(2003). Magnetic helicity is an important characteristic 
that strongly infiuences the PMF dynamics. Magnetic 
helicity is a conserved quantity during the subsequent 
evolution past the QCDPT. This leads to an inverse cas- 
cade producing magnetic fields at progressively larger 
scales. For this to work, it is important to know the 
magnetic helicity that is produced by the QCDPT. 

The magnetic helicity is defined as Jd^xA • 
with .B = V X A . In the domain wall model of 
IForbes fc Zhitnitskvl (j2000D . the magnetic helicity den- 
sity Hm is 

nM = A-B = A,B, , (5) 

for a PMF in the z direction, as discussed above. Because 
of strong CP violation during the QCDPT, magnetic 
helicity is produced through the alignment of magnetic 
and electric dipole moments of the nucleons. Thus, the 
electric field satisfies Ez — Bz (see IForbes fc Zhitnitskvl 
120001 ). From Maxwell's equations in the Weyl gauge we 
have 



E 



IdA 

~^~dt 



or 



-EzT , 



(6) 



where t ~ 1/Aqcd is the timescale for the QCDPT. 
From Eqs. dS]) and ^ one finds 



n 



=i(0.22 X 10^^ G)V(0.15GeV) 



(7) 



where we have assumed statistical homogeneity, so the 
result is gauge-independent. 

2.2. Comoving Values of Magnetic Field and Helicity 

The simple dilation due to the expansion of the Uni- 
verse significantly reduces the amplitude of both the MF 
and the magnetic hehcity created during the QCDPT. 
Defining and ao as the scale factors at the time of the 
QCDPT and today, respectively, we have 



ao 



0.15 GeV \ /15 



1/3 



(8) 



with g^=15, T^=0.15 GeV. 

The comoving (present) value of the PMF field Bin (the 
subscript "in" indicates that the QCD field is an initial 
PMF for further developed MHD dynamics) is given by 



(QCD) 



Si, 



which results in 
2 X 10"^ G. 



(9) 

As in the case of the PMF amplitude, magnetic helicity 
density experiences dilution due to the expansion of the 
Universe. The comoving (initial) value of the magnetic 
helicity density is given by 



M,in 



do/ 



10 



-39 



(G^ -Mpc). (10) 



This value is extremely small, and it is almost xlO^^ 
smaller than the maximal allowed magnetic helicity (see 
below). Such a small value of magnetic helicity density 
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is d ue to the thickne ss of the magnetic wah; see Eq. (|2]) 
and iKisshngerl ()2QQ3l ). On the other hand seve ral studies 
indica te strong CP violation during QCDPT (jKharzeevI 
[20M iVoloshinI [2003 : ICreutzl [20TTI ). In this case mag- 
netic helicity can reach its maximal value, if we assume 
the MF to be correlated over the Hubble scale \h^^ as 
will be explained below. Being more conservative we as- 
sume that the MF correlation should coincide with the 
bubble size ^m, see Sec. 3. The resulting magnetic he- 
licity will then be smaller than the maximal one by a 
factor of the order of ^m/A//^. As we will see in Sec. 4, 
the duration of the process is long enough to ensure that 
the maximal value of magnetic helicity is reached during 
the subsequent evolution. 

3. MAGNETIC FIELD SPECTRUM 



Follo wing our earlier studies (Kahniashvili et all 
[20103 ) we treat the initial PMF energy density £m as 
magnetic energy density injected into the cosmological 
plasma at the comoving length scale Aq which corre- 
sponds to the QCD bubble size. We recall that the 
PMF has been generated on the thin surfaces between 
colliding bubbles, while the correlation length scale of 
this PMF might be associated with the bubble length 
scale. In the following, we assume that the PMF spec- 
trum in Fourier space is sharply peaked at /cq = 27r/Ao. 
After generation, the PMF evolution (during the PT) 
depends sensitively on the length scale under consid- 
eration and on the presence of magnetic helicity (see 
"Harrison 1970; Turner & Widrow 1988; Cornwah 1997; 
Jackiw Pi 2000; Garrets on et al. 1992; Fi eld &: Carrolll 
20001 : [cTovannini 2000: iVachasp ati 2001J iCampanellil 



2007L 120091 : iDurrer et al.il201lL for magnetic hehcity gen- 
eration mechanisms). The expansion of the Universe 



leads to additional effects, in particular to a faster growth 
of the PMF correlation length. A distinctive effect is the 
different time behavior of the PMF decay. 

In the cosmological context most important is the 
difference between the growth of the comoving length 
scale (L ex a) and the Hubble radius {H~^ ex t, where 
t is physical time). This leads to additio nal effects 
in the PMF evolution (and damping) (s ee iSonl 119991 : 
iBaneriee & Jedamzi^ 120031 : ICaprini et all l2009f ) . Note 
that, to describe properly the dynamics of the pertur- 
bations in the expanding Universe, it is appropriate to 
switch to comoving quantities and to describe the pro- 
cesses in terms of conformal time 77 (iBrandenb urg et "aTl 
[T996). After this procedure the MHD equations include 
the effects of the expansion while retaining their conven- 
tional flat spacetime form. To keep the description as 
simple as possible we work with dimensionless quanti- 
ties, such as the normalized wave number^ 7 = Xo/H~^ 
and normalized energy density defined below. 

The coupling between the PMF and the plasma leads 
to a spreading of the fixed scale PMF over a wide range 
of length scales, thus forming the PMF spectrum. Af- 
ter a few turnover times the modified PMF spectrum is 
established (see Sec. 4 for details of the simulations). 

To show the coupling between the initial PT-generated 
MF and the plasma we give here the basic MHD equa- 

^ Here the subscript indicates again the moment of the PMF 
generation. 7 can be associated with the number of PMF bubbles 
within the Hubble radius, oc 7^. This value depends on the PT 
model: for the QCDPT we assume 7 ~ 0.15. 



tions for an incompressible conducting fluid (jBiskampI 
M) 



d_ 
df] 

d_ 
df] 



i; = (6- V)6- Vp + /k,(11) 
6=(6-V)i; + Vx/m, (12) 



with V • 6 = 0, where 77 is the conformal time, v(x^r]) 
is the fluid velocity, b{x^r]) = B {x , r]) / VAttw is the nor- 
malized MF, fxix^r]) and fui'^^Vj) are external forces 
driving the flow and the magnetic field {fx = /m = 
for the results presented below, but /m 7^ for pro- 
ducing initial conditions), v is the comoving viscosity of 
the fluid, A is the comoving resistivity, w = p ^ p \^ the 
enthalpy, p is the energy density, and p is the pressure 
of the plasma. Here we are interested in the radiation 
dominated epoch. 

To proceed we derive the Fourier transform of the PMF 
two point correlation function as 

{b*{k,rj)bj{k',v+r)) = {2nf5{k-k') F,^{k,T) ./[«(fc),T], 

(13) 

Such a presentation allows a direct analogy w ith hydro- 
dynamic turbulence CLand au fc Lifshitzl 119871 ). In fact, 
bi represents the Alfven velocity. The normalized energy 
density of the PMF is then £m = (&^)/2, while the ki- 
netic energy density is Sk = and the spectral 
correlation tensor is 



Kfik,r) 
(27r)3 

Here, Pij{k) 



+ iSijiki 



Sirk^ 



(14) 



ij — kikj /k'^ is the projection operator, 
5ij is the Kronecker delta, k = \k\, Siji is the totally an- 
tisymmetric tensor, and f<i{k) is an autocorrelation func- 
tion that determines the characteristic function f[f<i{k)^r] 
describing the temporal decorrelation of turbulent fluc- 
tuations. The function HMik^Vj) is the magnetic helic- 
ity spectrum. Note that EM{k) = /c^Pe(A:)/7r^, where 
PB{k) is the MF power spectrum. 

The power spectra of magnetic energy EM{k,r]) and 
magnetic helicity HM{k,r]) are related to magnetic en- 
ergy density and helicity density through Euij]) = 
dkEMikjT]) and Huiv) = ^kHM{k,r]), respec- 
tively. The magnetic correlation length. 



f 

Jo 



dkk ^EM{k,r]), 



(15) 



corresponds to the largest eddy length scale. All config- 
urations of the MF must satisfy the "realizability condi- 
tion" (Biskamp 2003 ) 



\nM{r])\ <'2^M{r])£M{r]). 



(16) 



Also, the velocity energy density spectrum Exik^r]) is 
related to the kinetic energy of the turbulent motions 
through fx (7/) = dkExik^T]). 

One of the main characteristics of the PMF is the cor- 
relation length and its growth. The maximal correlation 
length ^max for a causally generated PMF cannot exceed 
the Hubble radius^ at the time of generation H~^. The 

^ T he inflation generated PMF (Turner & Widrow 1988; R atral 
119921 ) correlation length can exceed the Hubble horizon today. 
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comoving length corresponding to the Hubble radius at 
generation is inversely proportional to the temperature 



Xh, = 5.3 X 10"^ Mpc 



0.15 GeV 



1/6 



(17) 



and is equal to 0.5 pc for the QCDPT with = 15 and 
= 0.15 GeV. 

The PMF spectrum is characterized not only by its 
spatial distribution, but also by its characteristic times: 
(i) the largest-size eddy turnover time tq — Iq/ va (where 
va is the r.m.s. Alfven velocity), which can also be used 
to determine the minimal duration of the source needed 
to justify the use of the stationary turbulence ap prox- 
imation (iProudmanI 119521 : iMonin fc YadomI I1975D : (ii) 
the direct cascade time-scale of the turbulence Tdc; and 
(iii) the large-scale turbulence decay time rig. 

The temporal characteristics of the MHD turbulence 
is given through the form of r), which is due 

to t he complex process o f MHD turbulence decorrela- 
tion (|Terrv fc Smith|[2"007[ ) and is currently not fully un- 
derstood. To pro ceed we employ Kraichnan's approach 
(|Kr aichnani 1 1964) and specify the decorrelation function 
/dc[/^(A:ph), t] = exp[— 7r/T:^(/cph)T^/4] defined within the 
inertial range, < k < kd- Here r is the duration of the 
turbulence process and /^(/cph) = ^M^^ph^/A/^Tr, where 
/Cph is the physical wavenumber related to the comoving 
k through /cph(a) = kao/a (ao is the value of the scale 
factor now), and em is the proper dissipation rate per 
unit enthalpy. Hence, we have (Kahniashvili et al. 2011) 



fdc[k,r] = exp 



2^ 
9 



k^/' 



(18) 



Here, k = k/ko is the normalized wavenumber and tq 
corresponds to the largest eddy turnover time. It is clear 
that after switching off the forcing, the turbulent motions 
are decorrelated within a few turnover times, and are in 
fact irrelevant to infiuence the large-scale MF. 

4. GROWTH OF CORRELATION LENGTH IN HELICAL 
TURBULENCE 

To assess the importance of a small initial magnetic he- 
licity, we perform direct numerical simulations of decay- 
ing MHD turbulence with an initial MF of finite relative 
magnetic helicity using different values, and a correlation 
length that is small compared with the scale of the 
domain Ai. 

4.1. Simulation technique 

We solve the compressible equations with the pressure 
given hy p = pc^ , where = 1 / a/3 is the sound speed 
for an ultra-rel ativistic gas. Following our earlier work 
(jKahniashvili et al.l 12010a) . we solve the equations gov- 
erning equations for the logarithmic density In p, the ve- 
locity and the magnetic vector potential A, in the 
form 



Dlnp 
Dt^ 
Di; 



-V • V, 



— = JxB-cjVln/> + /v 
Drj 



(19) 
(20) 



dA „ 

— =vxB- 
or] 



TM 



(21) 



where D/Dt^ = d/dr] -\- v ■ ^ is the advective deriva- 
tive, /vise = ^ (V^v + • v + G) is the viscous force 
in the compressible case with constant u and Gi = 



SijVjlnp as well as Sij = ^{vij + vj^i) — ^SijVk^k be- 
ing the trace-free rate of strain tensor. Furthermore, 
J = V X B/AiT is the normalized current density. We 
emphasize that /m = 0, except for producing initial con- 
ditions, as explained below. 

The bulk motions are always slow enough, so compress- 
ibility effects are not important. Similar to before, we ex- 
press the magnetic field in Alfven units, but now based on 
the volume average enthalpy, i.e., h = B/ ^At:{w)^ where 
w = I p for an ultrarelativistic gas. We use 512^ mesh- 
points in a domain of size (27r)^, so the lowest wavenum- 
ber in the domain is /ci = 1. We choose v = r] = in 
units of Cs/ki. 

4.2. Initial conditions 

A suitable initial condition is produced by simulat- 
ing for a short time interval (At ^ O.SAi/c^) with 
a random ^-correlated magnetic forcing term /m in 
the evolution equation for the magnetic vector poten- 
tial. This forcing term consists of plane monochro- 
matic waves with wavenumber k^ and fractional helicity 
(/m • V X fM)/{koflt) = 2a/{l + a^); in the following 
we quote the value of a. This procedure has the advan- 
tage that the magnetic and velocity fields used then for 
the subsequent decay calculations are obtained from a 
self-consistent solution to the MHD equations. 

4.3. Growth of helical structures 

In Figure [T] we show spectra of magnetic and kinetic 
energy, as well as the magnetic helicity scaled by A:/2, 
for a run with a = 0.03. Initially, kHM{k,r])/2 is well 
below the value of EM{k^ rj). However, at later times the 
two approach each other at large scales. This shows that 
the relative magnetic helicity increases during the decay. 
For the four times shown in Figure [TJ the rms Mach num- 
ber, Vrms/Cs^ is 0.05, 0.025, 0.012, and 0.007; 6rms/'^rms 

is around 3.4, and the Reynolds numbers Vrms£,M/^ are 
roughly 270 all cases. 

The growth of turbulent structures is particularly clear 
in the magnetic field (Figure [2]). The magnetic field 
drives correspondingly larger scale structures also in the 
velocity field. However, there are also strong small-scale 
fiuctuations in the velocity field that are not visible in 
the magnetic field; see the second row of Figure [2l 

In agreement with earlier simulations, we find that 
at small scales the magnetic energy is re-distributed 
by a direct cascade with a Kolmogorov-type spectrum, 
EM{k) (X k~^^^. At large scales a Batchelor spectrum^, 
Euik) (X /c^, is es tablished, which was used as initial con- 
dition already in Brandenburg et al. (1996). This spec- 
trum is in agreement w ith the analytical description of 
iDurrer fc Caprinil ()2003l ) who derived this result from the 
requirement of causality and the divergence free condi- 
tion. The earlier study of lHoganI ([1983), which thus vio- 
lates causality for magnetic energy, yielded a white noise 

Sometimes this spectral distribution, £^m(^) oc fc^, is called a 
von Karman spectrum (,PoDe..20Q0. ). 
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Fig. 1. — Spectra of magnetic energy (solid lines), kinetic energy 
(dotted), and magnetic helicity scaled by /c/2 (dashed) for a run 
with cr = 0.03 at three different times. At early times, HM(k,r]) 
can be negative at small values of which explains why the dashed 
line terminates in those cases. 

spectrum E{k) oc k'^ (Saffman spectrum) which we do 
observe for the spectral distribution of the kinetic energy 
EK{k) (X P. 

AA. Growth of turbulent length scales 

The evolution of magnetic correlation length and mag- 
netic energy during the MHD turbulence decay can be 
described using two indices and ue'- 

. .. 



(22) 



^m(^) Snivo) — 



(23) 



In this case, we can model the spectral energy density 
of the PMF using time-dependent large- and small-scale 
ranges: 



k^ when k < ki{r]) .^4) 
k~^^^ when k > ki{r]) ' ^ ^ 



where k = k/kj and ki{r]) = 27r/^M(^)- Hence, the evo- 
lution of the spectral amplitude £^0 foi" a given magnetic 
field spectrum will be [see Equations (|22]) and (|23|) ]: 



lYTT \r]o 



(25) 



Magnetic helicity crucially affects the e volut i on of 
the PMF (^Biskamp Miillerl [Tool [200(1 ISonl [Tool 



iChr istensson et alJ '20011, 120051 : iBanerjee fc JedamzikI 
[2003, 2004; Campanelh 200^. If the PMF has been gen- 
erated with small magnetic helicity, there are two main 
stages during the development of the MF spectrum: dur- 
ing the first stage (sometimes called direct cascade) the 
PMF dynamics is similar to that of the non-helical MF. 
The energy cascades from large to small scales where it 
decorrelates and dissipates: this is a standard forward 
cascade development. Since magnetic helicity is con- 
served, its fractional value increases and thus the end 
of this first stage is characterized by releasing turbulence 
to a maximally helical sta te (jBanerjee fc Jedamzikll2003l : 
IChristensson et aJ] l2005f ) when the realizability condi- 
tion (p!6]) is reached, the inverse-cascade stage starts. 
The conservation of magnetic helicity implies that the 
magnetic energy density decays in inverse proportion to 
the correlation length growth during the inverse cascade. 
The realizability condition implies that 



(26) 



so there is a minimum value for the correlation length. 
In Figure [3] we plot ^m{v) and Q^"(7?) for a = 1, 0.1, 
and 0.03. It turns out that, especially in the latter case 
with a = 0.03, the increase of remains slow (~ rj^^'^) 
as long as ^m{v) ^ ^m"(^)- However, since 1~Lm is es- 
sentially constant and Sm decreases approximately like 
77"^, the value of ^]S^"(^) soon reaches ^m(^)- When that 
happens, the field is essentially fully helical and the cor- 
relation length and the magnetic energy density evolve 
according to ^ v'^^^ and £m ^ 7^"^^^, respectively. 
Hence, we recover two distinctive phases in the MHD 
turbulence decay process: evolution of a weakly helical 
turbulence with = 1/2 and ue = —1, and fully he- 
lical turbulence with = 2/3 and ue = —2/3. Note 
that in the latter case Eo{r]) oc ^m(^)^m(^) = const 
[see Eq. ([25]) ] and the inverse cascade develops. Our re- 
sults ar e in excellent agr eement with Biskamp & Mi illerl 
(11999^, 'BiskamD & Miiller (2000), Banerjee & JedamiiS 
(2003), and Campanelh (2007). The dynamical process 
of PMF coupling with the cosmic plasma stops at the 
moment of r ecombination after which the PMF develops 
more slowly ([Brandenburg et "all 1 199 71 ). 

To calculate the time r^fuiiy when a fully helical state is 
reached, we only need to know the initial values ^m(^o) 
and ^m"(^o)- Since the latter approaches the former like 
7/^/^ the result is r/fuiiy = Vo[^m{vo)/^m'''{vo)? • Thus, in 
terms of the initial values of £m and Hm, a fully helical 
state is reached at the time 



^fully 



(27) 



Note that this time increases quadratically with decreas- 
ing initial value of Hm- In case of the strong CP violation 
during the QCDPT, when the initial magnetic helicity 
can reach values that are only £,m/^h^ times less than 
the maximal one (see Sec. 2), we get r^fuiiy = ^o/t^- 

4.5. Observed Magnetic Fields 

Galactic and cluster MF s are usua lly measured 
through Faraday rotation (see IValli3l2004l ) and, as men- 
tioned above, the value of the coherent magnetic field is 
of the order of a few /iC with a typical coherence scale of 
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Fig. 2. — Visualizations of Bx (upper row) and Vx (lower row) at three times during the magnetic decay of a weakly helical field with 
r = 0.03. 
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Fig. 3. Evolution of ^M(r]) (solid) and Cm^'C^) (dashed) for 
a = l (black), 0.1 (blue) and 0.03 (red). 

lOkpc,^ and cluster MFs have lower limits of the order 
of 10~ ^ G, and at least a few nG, with similar coherence 
scales (Clarke et al. 2001) and additional lower limits on 
the steepness of the magnetic power spectrum in clus- 
ters. Furthermore, simulations starting from a constant 
comoving magnetic fields of 10~^^ G suggest that mag- 
netic field generation in clusters can be sufficiently strong 
to explain Faraday rotation measurements (fPolag et al.l 
[2002; Jedamzik & Sigl 2011). 

^ Strong MFs have been detected through Faraday rotation of 
distant quasars proving that the MFs comparable to those observed 




Fig. 4. — Spectral energy density of the turbulent magnetic field 
\ogEM(k) (color coded) in a representation of magnetic correla- 
tion length versus temperature. The thick solid line shows the 
evolution of the magnetic correlation length ^m(T). The magnetic 
correlation length starts to grow after the QCD phase transition at 
Ti, = 0.15 GeV, when = 0.075 pc. The transparent dashed area 
corresponds to decorrelated magnetic field. White arrows show the 
direction of the evolution during the expansion of the Universe. 
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-1/2 and ue = 1- 



today are seen at high redshift z ^ 3 



. .2008.) . 



Figure |4] shows the spe ctral energy de nsity of the 
QCDPT-generated MF (see lKisslingerll20Q3f ) with respect 
to temperature and correlation length in weakly helical 
turbulence. Initially the integral scale of the MHD tur- 
bulence is set by the QCDPT bubble scale (lower right 
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Fig. 5. — Similar to Fig.0] but for the case in which the magnetic 
field reaches a fully helical state within the considered expansion 
time of the Universe. Initially, during the growth of magnetic he- 
licity correlation length (solid line) and energy evolve according 
to = —1/2 and ue = 1- After reaching a fully helical state, 
correlation length (dashed line) and energy evolve according to 
= —2/3 and = 2/3. 

corner of the diagram). The thick sohd hne marks the di- 
vision between the evolution of large-scale (plain colored 
region) and small-scale (hashed colored region) magnetic 
fields. White arrows indicate the direction of the evo- 
lution during the Universe expansion. At scales below 
the integral scale of the turbulence, the magnetic field 
undergoes exponential decorrelation; see Eq. ([18]). The 
integral scale of the MHD turbulence increases, reaching 
= Ikpc at T = leV. Here we have used = 1/2, 
He = —1, with initial magnetic helicity corresponding to 
that given by Eq. (p!Q|) . 

Figure [5] shows the spectral energy density of a 
QCDPT-generated MF in the case when the initial heli- 
cal turbulence reaches the fully helical case during the 
expansion of the Universe with 7^fuiiy/^o = 1/(0.15)^. 
The thick solid line marks the evolution of the magnetic 
correlation length until the magnetic helicity reaches its 
maximally allowed value. In this time interval the de- 
cay law for weakly helical turbulence with = 1/2 and 
He = —1 is applied. After the time r^fuiiy when maxi- 
mal magnetic helicity is reached, the correlation length 
follows the black dashed line and the MF evolution fol- 
lows that of the fully helical case with = 2/3 and 
He = —2/3. The integral scale of the MHD turbulence 
reaches = lOkpc at T = 1 eV. 

The presented model is somewhat idealized since it 
ignores the time of Silk damping due to large corre- 
lation lengths for pho ton and neutrino viscosity (see 
iJedamzik fc Sigll 1201 1[ ). This is justified since it only 
delays the evolution but does not destroy the field 
(jBrandenburg et aTl 119971 ). Therefore, we can present 
here only upper values for QCDPT MFs within the model 
bylKisshnger (2003). 

The final amplitude of the MF can be estimated 
through two different approaches, (i) We compute the 
total magnetic energy density, i.e. £m = EM{k) 
and make the assumption that all energy is again given 
only at one scale that corresponds to the integral scale 



at this moment, i.e. ^eff = V SttSm • (ii) Another ap- 
proach is to compute the strength of t he MF , B(X\ at 
a given scale A. Since observations (|Valle e' "2004") do 
not allow us to properly reconstruct the configuration 
of the MF we adopt first a n "effective" MF approach 
(see Kah niashvili et aTll2011[ ). The resulting value of the 
effective MF in our model of weakly helical turbulence 
with = 1 kpc reaches 5 x 10~^ nG, while in the case 
of a fully helical configuration with = 10 kpc we find 
7 X 10-^ nG. 

5. CONCLUSIONS 

In this paper we have considered QCDPT-generated 
PMFs and their evolution in an expanding Universe ac- 
counting for the effects of MHD turbulence to explain the 
seed MFs of clusters and galax ies. We consider the MF 
generation model proposed bv iKisslingerl (|2003f ). which 
yields an initial state of weakly helical MHD turbulence. 
We also study the possibility of strong CP violation ac- 
cording to jForbes fc Zh itnitskv (2000), which yields an 
initial state with much higher magnetic helicity at a time 
when maximal helicity of the MHD turbulence is reached 
during the expansion of the Universe. The initial seed 
MF is generated via QCDPT bubble collisions with a co- 
moving correlation length of the order of 0.1 pc and with 
a comoving amplitude of the order of 20 nG. The initial 
magnetic helicity is determined by the thickness of the 
surface between two colliding bubbles and is extremely 
small if no strongs; CP violation i s assumed (|Kharzeevl 
[2QQ6I : iVoloshinI [2QM : 'Creutz" 'MT) . During the expan- 
sion of the Universe there are different processes that 
affect the correlation length and the strength of the MF: 
first of all, during the PT the field is initially peaked at a 
given scale and then spreads out within a wide range of 
wavelengths, establishing a Kolmogorov-like spectrum, 
EM{k) oc A:~^/^, at small scales and a Batchelor spec- 
trum, EM{k) oc /c^, at large scales. If the PMF was gen- 
erated without being maximally helical, the magnetic he- 
licity experiences a steady growth. One of the results ob- 
tained in this paper is an estimate of the timescale within 
which the field starts to be fully hel ical. In the case of 
an extremely weakly helical field (Kisslingerl 120031) . the 
available time to produce a fully helical PMF may be too 
long. The growth of the correlation length follows then 
a oc T~^/^ law. For moderate or reasonably small 
initial magnetic helicity (even for a > 10~^ — 10~^), the 
evolution timescale is long enough so that during the first 
stage of evolution, magnetic helicity grows to its maximal 
value. During the next stages (after magnetic helicity has 
reached its maximal value) the correlation length experi- 
ences a steady growth with the scaling law oc T~^/^ 
while the energy density is decreasing in the opposite 
way keeping magnetic helicity almost constant. Finally, 
at recombination the growth of the correlation length 
slows down. The resulting correlation length in the most 
optimistic scenarios is around 10 kpc and the amplitude 
of the MF is around 0.007 nG. Assuming t hat the MF 
is ani plified during the growth of structures ( Dolag et al.l 
I2OO2I ) , such a field might well be strong enough to explain 
the observed MF in galaxies and clusters. On the other 
hand, observations of the CMB fiuctuatio ns are sensitive 
to PMFs of the or der of a few nG (see I Shaw fc LewisI 
I2OIOI : lYamazaki et al...2010. . and references therein). 
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Another possible signature of QCDPT-generated 
magnetic fields i s a gravitational wave signal 
(iKahniashvili et aP l2QlQbl ) that might be indirectl y 
detected through pulsar timing (jPurrer et al.l l2Qllf ). 
The gravitational waves signal from PTs is usually 
computed assuming short duration of the source (either 
turbulence or PMF anisotropic stress). On the other 
hand, due to the free decay of MHD turbulence, the 
source of gravitational waves acts also after the end of 
PTs. For short duration sources, the peak frequency of 
the gravitational waves is fully determined by the source 
characteristics. In particular, for QCDPT-generated 
gravitational waves it is far too weak to be detected 
though gravitational waves via ground or space based 
missions. Long duration sources might in principle 
substantially change the peak frequency as well as the 
amplitude of the signal. We plan to address this issue in 
future work. 
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